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ABSTRACT 


The paper outlines the theoretical basis behind the defini- 
tion of frequency stability in the time domain. Various 
types of variances are examined. Their differences and 
interrelation are pointed out. Systems that are generally 
used in the measurement of these variances are described. 


INTRODUCTION 

Radio frequency sources give an output signal which, in general, 
is affected by small fluctuations in their amplitude, phase or frequency. 
The nature of these fluctuations may be random or deterministic. Due to 
the large number of users of these sources and the variety of fields in 
which they are applied, a problem has arisen in the method of character- 
ization of the frequency fluctuations. In the case of the deterministic 
fluctuations such as linear frequency drift, it has been found, in most 
cases, that a specification of the fractional frequency deviation on a 
per day or per month basis is satisfactory. However, in the case of non 
deterministic or random fluctuations, a mathematical treatment based on 
probability concepts is necessary. In the past, depending on the field 
of interest, several methods of characterization have been used and some- 
times have led to confusion. 

More recently a proposal has been made by the IEEE Subcommittee 
on Frequency Stability that the spectral density S (f) of the fractional 
frequency fluctuation, y, and the two sample variance, a y (t), be used to 
characterize frequency stability respectively in the frequency domain and 
in the time domain [l] . These two parameters have resulted as a logical 
conclusion from a rather large amount of theoretical and experimental 
work on the subject [2], [3]. However, other parameters, specially ir. *-he 
time domain, have been studied and have been found most interesting in 
specific cases [ 4 ] , [ 5 ] , [6] , [ 7 ] . The subject still raises great interest. 


247 


In the present paper we give a brief description of the main theo- 
retical concepts involved in the definition of the frequency stability in 
the time domain. The interrelation through the spectral density S y (f) 
of some of the various variances that have been studied up to now is made 
explicit. Finally several systems used for the measurement of stability 
in the time domain are described. 

This paper was prepared at the request of the Program Committee of 
this Meeting. Several excellent reviews and tutorial papers on this sub- 
ject have been published in recent years [8] , [9] , [lO], [ll] , [l2] . Conse- 
quently it appears difficult in the writing of such a tutorial introduc- 
tion, to avoid repetitions or to improve on all these papers. The most 
one can do at this stage, is to present the material in a slightly diffe- 
rent manner. In particular the present text, specially on the theoreti- 
cal section, owes a great deal to he recent review by Dr. Rutman [8]; 
the reader is strongly encouraged to consult that excellent article. 

A - THEORY 


Definitions 


According to a well accepted notation, the instantaneous output 
voltage of a signal generator can be written as [l]: 

V(t) = [v 0 + e(t>] sin [ 2tt v„t+«(t)] , (1) 

where Vo and Vo are the nominal amplitude and frequency and e(t) and 
?(t) are amplitude and phase fluctuations. It is assumed that e(t) is 
very small with respect to Vq and can be entirely neglected. The ins- 
tantaneous frequency of the oscillator is defined as: 


v(t) = v 0 ( 

where 'Q(t) stands 
quency fluctuation 


1 + 


*(t) i 
2irvo ' 


t 


for dif(t) / dt . 
as: 


( 2 ) 


We also define the fractional fre- 


y ( t ) = till 

’ 2ttv 0 


and we assume that 

« 1 . 


till 

2itvo 


(3) 


(4) 


A proposed means to characterize the frequency stability of an 
oscillator is the spectral density of y denoted by Sy(f) . Its dimen- 
sions are Hz 1 . A measurement giving an estimate of Sy(f) would then 
characterize the stability of the oscillator in the frequency domain. 
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This can be done in practice with a spectrum analyzer. However, one must 
be aware that only estimates of S (f) can be obtained because of fre- 
quency range limitations and finite observation times. 

Experimental studies of various frequency sources have shown that 
for all practical purposes, the frequency fluctuations spectrum of the 
most common oscillators can be represented by a truncated polynomial in 
the Fourier frequency domain: 

S y (f)-h“f“, (5) 

where a is an integer, ranging from -2 to +2. The frequency fluctua- 
tions spectral density can be related to the phase fluctuations spectral 
density through the relation 

V £) - V f > • < 6 > 

The a has been associated with various types of fluctuations 
either in the phase or frequency representation. 


f 


< 




a 


Type of fluctuations 


2 

1 

0 

-1 

-2 


White Phase 
Flicker of Phase 
White Frequency 
Flicker of Frequency 
Random Walk of Frequency 


However, another characterization of the frequency stability can be made 
by considering that y(t) is a random function of time. Then, a statis- 
tical parameter measuring in some sense the excursions of the values y 
of the random function y(t) around its mean value should characterize 
the frequency stability of the oscillator. In statistics, the standard 
deviation 0 or the variance 0 2 Is often used as statistical parameter. 

We could define a variance for the instantaneous frequency. In practice 
the frequency is measured over a time interval ,t, called the averaging 
time and the variance is calculated through the relation: 


0 2 (t) 



(7) 


where < > means an average over an infinite number of samples. Without 
going too deeply into questions of sta. .stical concern, we make a few 
assumptions. First, we assume that stationnarity applies to our model. 
By this we mean that a displacement of the time coordinates does not 
change the statistics of our ensemble. Secondly we assume ergodi- 
city, that is, averages over the ensemble can be replaced by time aver- 
ages on one of the samples. The sign < > in equation (7) then becomes 
a time average. Furthermore, we can assume that <y> « 0 and o 2 (x) 


i 


ii 
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becomes simply < ( y ) > . Since < ( y ) > Implies an infinite time aver- 
age this variance is an ? ideal theoretical concept which is commonly cal- 
led the true variance I (t) . In practice it is clear that one can only 
do measurements either for a finite time or on a finite number of samples 
and thus obtain an estimate of this ideal statistical parameter. Further- 
more it is found that I z (x) diverges for certain tyoes of noise such as 
flicker frequency fluctuations. In order to avoid these problems, various 
scientists have proposed several types of variances obtained from limited 
amount of samples. As will be explained below, one of them, the "two - 
adjacent - sample variance" studied by Barnes and Allan, has been propo- 
sed as a time domain measurement of frequency stability [l] . 

An attempt at measuring time domain frequency stability 

Before defining the two sample variance let us examine the process 
of frequency measurement itself. We assume that a digital frequency 
counter is used to measure the frequency. The measurement is then made 
over a finite time x and one obtains an average of the frequency ove 
this time interval x . In other words, the counter gives the number of 
cycles during the time interval x . Figure 1 is an experimental 
arrangement by which the frequency of oscillator (1) is measured, oscill- 
ator (2) being used as the time base of the counter. For the purpose of 
simplifying the picture let us assume that oscillator (2) is perfect in 
the sense that its frequency is free of fluctuations. All fluctuations 
in the measurements would then come from frequency fluctuations of os- 
cillator (1) . The counter takes measurements in the sequence shown in 
figure 2(a). The result for N measurements may be as shown in figure 


Here x ■ t. . , - t, and 
k+l k 


V " v « < i +y kT > • 

(8) 



1 N 

<V>--r: )" V. 

N “ ix 

(9) 


The average value of the random variable v is only an estimate of the 
actual average frequency, average being done on N samples. One may then 
calculate for the N samples the variance: 

, 1 N 2 

- nTT l ( - < v > ) . (10) 

v j-i J 

In order to continue the analysis it is assumed that the oscilla- 
tor does not show syst* aatic drifts with time. If such drifts are pre- 
sent, they are. removed from the data and the following analysis applies. 

The variance of the random variable y may be readily written as: 
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1. First we note that a y (N, T, t ) is itself a random variable. 
It is an estimate of the true variance I 2 (t) made on N samples, its 
average value made c.. several sets on N samples, cr 2 (N, T, t ) , should 
be close to the true variance. At the limit whete T4 tends to infinity 
it should be equal to I 2 (t) . 


2 , The results of the experiments discussed above can be used to 
form an his gram as illustrated in figure 3. If the number N cf sam- 
ples is large enough the figure may be rather smooth and a good estimate 
of o 2 ( N, T, t ) can be obtained from this curve through a measurement 
of tne half width at half the height. (For a normal di- ; . ibution 

o - 1.17 ( ). 

* 

3. If the experiment is repeated for other averaging times t , 
the results obtained with the histogram technique may behave as shown in 
figure 4. The value of / o y 2 ( N, T, t ) for each of the histograms may 
then be plotted as a function of t . The results are shown in figure 5. 
These results appear very interesting, and give ■*- a sense an * idication 
of the frequency stability of the oscillator in the time domain. However, 
several difficulties arise when the technique of measurements is changed. 
For example it is observed that in the region where the variance varies 

as 1 /t 2 an increase in the number of samples N does not alter Its value, 
providing this number of samples is large enough. This is not the case 
in the region where a y 2 (N, T ) is independent of r . In that region 
an increase in the number of. samples N shows up by an increase of the 
variance. One is then faced with a problem of a variance whose va.urs 
depend on the number of samples. Furthermore, in the region whet 1 o 2 
varies as 1/ T 2 the value of o 2 depends on the frequency bandwidth of 
the measurement system. 


The two-sample variance 

In order to avoid these problems, ati„ to facilitate intercompari- 
son between the results reported by workers in various fields, one has 
tnen to make a choice on N and preferably the ratio T/t . It has been 
proposed, folltwing the work of Allan [l3J: 

1) that the following weighed sample variance be used: 

N N 2 

- nit l (y± l <V) S < 12 > 

i-1 j-1 J 


2) that N - 2 -* two-sample variance ; 

T 

3) that — =1 , or no dead time between measurements [l], T!--s 

variance is a random variable and its average, abreviato'* 
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(13) 


<jy 2 (x) , is given by: 

2 

o y 2 (x) = < o y 2 ( 2, T, t ) > = J < ( y k+1 - y k ) 

It is proposed as a characterization of frequency stability in the time 
domain . 



This variance has veij interesting properties: 

1. It is standardized in the sense that N and T are fixed. 

2. It is equal to the true variance I 2 (x) for white frequency 
noise. It is equal to 3/2 I 2 (x) for white phase noise thus close to it. 

3. It is convergent for all types of frequency noise normally 
encountered in oscillators, that is the five power laws mentioned earlier. 

4. Although by definition, one is still faced with an average on 
an infinite number of samples (in this sense a y 2 (x) is still an ideal- 
ization), good estimates of it can be obtained by a relatively limited 
number of measurements, m , of the pairs: 




y k+l - y k 


2 

) 

j 


(14) 


For n > 10 it has been shown chat the confidence interval to be 
associated to a (x) in such an estimate if of the order of [l4] . 


Confidence interval = K a (x) //m 

ay 


where Kq is a constant depending on the power law predominent, but is 
not far from unity. 


5. Finally, tables have been compiled which translates from one 
type of variance to another in relation to variations in N and the ratio 
T/t , and this, for the five power laws commonly encountered in oscilla- 
tors [l5j . Of particular interest is the bias function 


B 2 ( } , u ) 


<0 y 2 (2, T, x ) > 

< ( 2, x, x ) > * 


since in general a simple counter will be characterized by a dead time 
(T/x ^ 1) between successive measurements. 


The main disadvantages of a y 2 (x) are: 

1) it diverges for power law spectral densities greater than -2; 

2) It does not discriminate between white phase noise and flicker 
of phase noise. 
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Relation between the time domain and the frequency domain; other types 
of variance 


The time domain frequency stability, as characterized earlier 
either through the true variance or the two-sample variance, can be inter- 
preted in a different way. The operation of the counter, averaging the 
frequency for a time t may be thought of as a filtering operation. The 
transfer function , H(f) , of this equivalent filter is then the Fourier 
transform of the impulse response h(t) . It can be shown that the time 
domain frequency stability is then given by [16 ] , [l7] . 


a 2 ( N, T, T ) > = I S (f) | H(f ) | 2 df , 

Jo y 


(15) 


where S (f) is the spectral density of frequency fluctuations. In the 
case of^the true variance and of the Allan variance, we have: 


f 


Hj(f ) = 


1 

T 


. U)X 
sin ~ 2 ~ 

OJT 
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t < - T 


0 

t < -T 



1 

-T < t < 0 

-T < t < 0 , 

/—* \ 
Pt 
s-/ 

II 

7Tt 

i 

t>0 



0 < t < X 



VO 

t > X 



2d) T 



(16) 


sm 


H .(f) = 

a aa 2 


(17) 


These relations are illustrated in figure 6. 


This "transfer function approach" has been exploited by several 
authors to elaborate new types of variances for characterizing oscillator 
frequency stability in the time domain. It is the equivalent of digital 
filtering used in data processing. 

Hadamard Variance 


The sequence of measurements shown in figure 6 for the Allan 
variance, which consists of two samples (N=2) can be expanded to a se- 
quence of a greater number of samples . A sequence for the case where 
N = 10 is shown in figure 6(e), where the impulse response of the equiva- 
lent filter is plotted as a function of t . The variance for this se- 
quence is then [l8|: 

<o R 2 (n, t d , x) > - < (y 1 -y 2 + y 3 - ••• -y N ) > 2 > . ( 18 > 

where T R is the dead time between measurements. The square modulus of the 
transfer function of this filter is [20 J : 
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(19) 


V f » I 2 


f sin tt f > 2 r sin N tt T f 2 

*■ IT T f COS Ilf J 


It is illustrated in figure 6(f) for no dead time and N = 10 . The charac- 
teristics of this variance, interesting for the topic of frequency stab- 
ility characterization are as follows: 

a) For the case of no dead time (1^ = 0), the transfer function of 
the equivalent filter h'j a main lobe, centered at f^ = l/2t, and whose 
width is equal to: 

IT 2 

f = -Jr- — — (equivalent rectangular filter) 


It may thus be made very narrow by increasing N . In this sense this 
filtering process is well suited for spectral analysis and this property 
has been exploited to obtain the spectral density if frequency fluctua- 
tions using time domain measurements [20] . This is easily seen from re- 
lation (15) by realizing that | H(f) | 2 can be approximated by a narrow 
square window over which S (f) does not vary much. The spectral density 
is then given by: ^ 

S y (f l ) = N <0 H 2 <N ’ T * T } > (20) 

This is seen to lend itself to a straightforward computation, in order to 
obtain an estimate of the spectral density without the use of a spectrum 
analyser . 

b) However, one should be aware of severe limitations in this tech- 
nique. There exist secondary side lobes in H^(f ) , which appear at harmo- 
nics of fi=l/2x, for no dead time between measurements. These can be 
minimized t>y proper adjustment of the dead time between measurements or 
by proper weighing of the samples of the measurements. The properties of 
the transfer function have been well studied in the case of a weighing by 
the binomial coefficients and by a pseudo-sinusoidal function [l8j , [19] , 
[ 20 ]. ' ' ' ' 

Modified sample variance 

Boileau and Picinbono have introduced a variance which can be in- 
terpreted in the case of its digital realization as a modified sample 
variance. [6j Their relations, when translated into the notation generally 
adopted in the field of frequency stability gives a variance as follows: 

a mod ( N * T * T * =7 (N+1)/2~N * (21) 

where N is odd. Thus the central sample of the set of N 



samples. Rutman has examined the case when N “ 3 and T * x [8] . The mea- 
surement sequence and the square modulus of the transfer function are 
shown in figure 6(g) and 6(h) respectively. The variance is given by: 


< o 


mod 


( 3, x, t ) 


f 

■in 


„ 16 sin 6 it x f 

V } 9 (ttx f)2 


df 


( 22 ) 


The modified sample variance has the advantage of being convergent for all 
five power law spectral densities examined up to now, plus two others, 
where a = -3 and at = -4 . However, it does not discriminate between white 
phase noise and flicker phase noise better than o^(x) . 

A special case of the sample variance 

It is general practice to study the frequency stability or measure 
the variance as a function of averaging time x . This is what is done in 
the sample variance described above. In opposition to this practice, 

De Prins and Cornelissen [7] have studied frequency fluctuations over in- 
tervals T for fixed averaging times x . One can then in principle study 
the variance o 2 (T) of the instantaneous frequency (x-*0) as a function 
of the time interval T . This variance is very different from the one 
described previously. In the present case all values of y(t) are scan- 
ned as T is varied in opposition to the averaging over x considered up 
to now. 


In the case where x does not tend to zero but is fixed at a given 
value, the variance then becomes a special case of the sample variance 
and the behaviour of ct^ 2 (T) can be obtained from the bias function, B^. 

The high pass variance 

A close look at equation (15) suggests that a 2 (x) can actually be 
defined through the transfer function H(f) of the equivalent filter cor- 
responding to the measurement sequence. Rutman has suggested that this 
approach could be taken even if the actual measurement sequence was not 
existing [5]. Then, H(f) could be given the shape desired. Of course 
the inverse Fourier transform of H(f) is not necessarily a step wise 
function that could be implemented in a straightforward manner by a 
counting technique. Other measurement techniques have then to be imple- 
mented. 


In this approach the variance is written: 
f°° 

a2(T) = 72 J o M f ) I M f ) I 2 df . (23) 

where S^(f) is the phase spectral density and is related to S y (f) 
through relation (6). The variance is then defined in terms or H^(f) 
and not in terms of the measurement sequence. 
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From equation (23) one sees that for the Allan variance, the square 
modulus of the phase transfer function of the equivalent filter is 

| H^(f) | = sin 1 * it f x . (See figure *.7)) (24) 

This is essentially a high pass filter, having an oscillating nature with 
a period 1/x . Low frequency components f < (irx) -1 are filtered out, this 
being an essential character of the Allan variance. Consequently it ap- 
pears that | H^(f) | 2 could be essentially a high pass filter and essen- 
tially the same character for the variance would be obtained. In fact, 
calculations show that, when a second order high pass filter, with cut 
off frequency f c = (irx) -1 , is used to calculate a so called "high pass 
variance", the general behaviour with the power law spectral density is 
essentially the same as the behaviour of the Allan variance. Both vari- 
ances have the same asymptotic slopes with x and both variances cannot 
differentiate between white phase noise and flicker phase noise. 

Band pass variance 

Following this line of thought and recognizing the nature of the 
limitations of the high pass variance, Rutman [8] has suggested that a 
bandpass filter be used for | H^(f) [ 2 with a center frequency equal to 
(l/2x) and a constant Q factor, say equal to 1. In that case the beha- 
viour of a p (x) is quite different from that of o y (x) or ©^(x) ; it 
shows complete discrimination between the five power law models in its 
asymptotic behaviour as a function of x . 

Of course, the method for measuring ° B p( x ) is not a conventional 
one incorporating a frequency counter. One uses a phase comparator 
(loose-phase-lock technique), a bandpass filter and a r.m.s. voltmeter. 

In this sense, it is the same type of system as the one used in reference 
[l7] and essentially falls in the class of systems used for studying fre- 
quency stability in the frequency domain. It appears natural to think 
of frequency stability measurements in the time domain as being done 
through a time sequential technique and a statistical analysis of the 
resulting data. This should be kept in mind in the practical implementa- 
tion of systems designed for the measurement of frequency stability in 
the time domain. 

An unified approach 

In the previous sections, various types of variance were examined. 
The approach taken has been one in which the measurement sequence was 
identified; the transfer function of the equivalent filter implementing 
the impulse response for the sequence in question was established, and 
the variances could be calculated through relation (15) . This method is 
very useful in pointing out the limit of utilization of a particular va- 
riance in respect to the power law frequency model and also in understan- 
ding the reason of these limits. 
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These variances however have all been introduced as particular 
cases for special needs. Recently, Lindsey, Chie, Leavitt and Lewis [2l] , 
[22], [23] have introduced in the picture a method of analysis called the 
"structure function approach", which emphasizes the fundamental ties bet- 
ween these variances rather than their differences. 

The k 1 -* 1 average frequency fluctuation over time t can be written 
as a difference of phase: 

*k' 2^[* t k +,) • (25) 

We may define the first difference or first increment in phase as: 

A (1) ^(t k , t) = *<t k +T) - *(t k ) . (26) 

The frequency difference (7 -y, ) appears as a second difference in 

• K. * X 1C 

phase: 

^k+i -7 ^ ‘^7 r -'« t k +T) ' + <(t k- T) J • <27) 

and the second phase inc-ement is defined as: 

A (2) ^(t k , t) = ^(t k + x) - 2*(t k ) + ^(t k -T) . (28) 

It is readily realized that the first difference in phase is used in the 
definition of the true variance, while the second difference in phase is 
used in the definition of the Allan variance (two samples). Lesage and 
Audoin [20J have proposed to continue the process further and have obtai- 
ned an expression for the n fc h difference in phase, which includes the 
binomial coefficient as a weighing factor. This analysis has led them to 
the implementation of the Hadamard variance in which the measurement se- 
quence is weighed by the binomial coefficients. 

From this it appears that a common basis may be expected under the 
definition of the various variances examined above. The rank of the phase 
increment appears to play a major role. In the approach of Lindsey and 
Chie this point is stressed. The phase increment is defined as: 

N 

A N 'e(f,T)= l (-1) K fN]i^(t+ (N-k) T) , (29) 

k=0 V K ; 

( N ) N' 

k J = k , T (Binomial Coefficient) . (30) 

The structure function of phase is then defined as: 

D (N) (t) - < (A (N) *(t, t) > . ( 31 ) 
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Stationnarity of the difference is assumed in the wide sense. 
It is then shown that the variances defined earlier can all be expressed 
in terms of these structure functions: 


True variance 
I 2 (t) * 

All an variance 

2 


1 »S 1, <t) , 


(2tt vq t) 2 X 3 


V (T) ■ 2 ( 21 . D “ )(T) : 

Modified three-sample variance 


<0 y 2 mod (3 - T - T)> 


9T2VW D ” >(,) = 


Hadamard variance (weighed by the binomial coefficients) 

j 2 

H.B.C. 


< °U D r ( N, T, T ) > - T o — d4 N) (T) ; 


( 2tt v 0 


(32) 


(33) 


(34) 


(35) 


This approach thus clearly shows that the variances utilized up to now by 
various authors have a common basis, in occurence, a structure function of 
phase. On the other hand, this structure function is related to the spec- 
tral density through the relation: 

D^(t) = 2 2(N - 1) (2ttv 0 t) 2 S y (f) S ~ -j yy df , (36) 


which effectively, as stressed earlier, provides means for evaluating the 
spectral density S y (f) as filtered with a transfer function 


H sf (f > 


, 2N , 
sin 7T i f 

(ir x f) 2 


2(N-1) 


(37) 


through sequential measurements in the time domain. 

In the previous paragraphs a structure function of phase was in- 
troduced as a means for describing frequency stability in the time domain. 
A structure function of frequency, however, can also be introduced to 
describe frequency stability; it is written D (x) . Lindsey and Chie 
[22] have shown, by studying the mathematical ^differences between these 
structure functions, that the true variance is essentially a measure of 
phase instability while 'the Allan variance is a measure of frequency in- 
stability. This type of reasoning has led them to suggest that a func- 
tion of the product of the two types of variance could be a parameter by 
which frequency stability in the time domain could be characterized. 
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It should be pointed out that the characterization of the stabil- 
ity of an oscillator could in principle be made through tables of the 

structure functions dW(t) . The user could calculate from these tables ’* 

the type of variance that is best suited for his particular application. 

In a sense, these structure functions can be thought of, as characterizing 
completely the frequency stability of the oscillator in the time domain 
in the same sense as S^(f) , does it in the frequency domain. 

Long time frequency fluctuations 

In real oscillators it is possible to observe very long term fre- 
quency fluctuations, that is, very slow fluctuations which may appear over r 

periods of days, months and years. These may originate either from slow 
random fluctuations or from deterministic drifts in the behaviour of the 
oscillator. 

Slow random fluctuations 


The above analysis was limited to five power laws of the spectral 
density (-2 < a < +2) . Very slow frequency fluctuations predominate at 
very low frequencies and are thus represented by more negative slope power 
laws such as a = -3 or -4 . The greatest negative slope that the Allan 
variance can handle is o = -2 , a random walk of frequency type of noise; 
for more negative slopes it diverges. Since direct spectral analysis of 
these slow fluctuations is not experimentally feasible, it appears that 
the other types of variance mentioned earlier may be useful. Intact the 
modified three sample variance 0r y 2 !nO( j( 3, x, x ) converges for f - and 
f -4 types of noise with respective slope i 2 and t 3 . Consequently 
for very slow frequency fluctuations one may then have to use a variance 
different from the two sample variance in order to have meaningfull in- 
terpretation of time domain data. 

Deterministic drifts 


Systematic drifts are generally observed in oscillators. These 
drifts may be represented by a polynomial [9j . A model for frac- 
tional frequency drifts is: 

y(t) = d^t 1 + d^t 2 + ... d n t n . (38) 

For the first term, rep.' - iting a linear frequency drift, the Allan va- 
riance varies as t 2 . ii ^ time dependent for higher order drifts. The 
other variances mentioned auove consequently are found very useful in 
characterizing these higher order polynomial drifts. In particular the 
modified three sample variance varies as t 2 for quadratic frequency 
drifts while ti»e Hadamard variance weighed by binomial coefficients varies 
as t 3 for cubic frequency drifts. The behaviour of the asymptotic value 
of o(x) versus t for various power laws of spectral density and various 
orders of frequency drifts are summarized in table 1. It should be under- 
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stood that the structure function approach can also be applied to analyze 
these long term frequency fluctuations, although emphasis has not been 
placed on this point in this paper. 

In this section we have examined the properties of several types 
of variances that can be used for characterizing frequency stability of 
oscillators in the time domain. Differences between these variances 
have been shown; their interrelation through tho fractional frequency 
fluctuations spectral densicy has also been emphasized. In fact, it has 
been shown that these variances, in the case of random fluctuations, are 
essentially elegant means of representing the spectral density S y (f) 
through parameters which can be measured with simple systems implemented 
with a frequency counter and a calculator for doing statistical analysis. 
Such systems will now be described. 

B - TIME DOMAIN MEASUREMENT SYSTEMS 


Ideally, frequency stability measurements require a frequency refe- 
rence much more stable than the oscillator to be studied. A system real- 
izing this condition can be implemented easily for the measurement of the 
frequency stability of most common oscillators. However, for stable, 
state-of-the art, oscillators, it is necessary that the reference oscil- 
lator be at least as stable as the oscillator studied. Erequency stab- 
ility measurements in the time domain can be done by two different methods 
which are related to the detection of two different parameters: frequency 
and phase. 

One can determine the mean frequency over finite observation times 
and calculate, for a given number of samples, a certain variance as des- 
cribed in the previous section. Although these measurements involve the 
well developed technology of frequency or period counting, their use is 
limited by the fact that they give information only on the mean frequency 
and the frequency fluctuations. Furthermore, in many systems, the mea- 
surement samples are not adjacent in time which, in jome cases, alter the 
value of the statistical parameter calculated. 

Many experimental set up's have been considered during the last 
several years, as the technology evolved 'see for example reference [24], 
[25j and [26] . We shall limit our discussion to a few of them illustra- 
ting their main principles. 

The most simple set up is the "Direct Frequency Counting System" 
shown in figure 8. In this case the reference osc-Lllator is the frequency 
counter time base. It is suitable for the study of low performance oscil- 
lators. The output of the counter gives the mean frequency over a prese- 
lected time interval. Its recording on a trip chart, magnetic tape or in 
a digital memory allows one to use best estimate curve fitting methods, 
in order to find any systematic frequency trends or drifts. These drifts 
are removed prior to statistical analysis. The two sample variance (Allan 
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variance) or any other desired variances are calculated from the set of 
corrected data. The main limitation of this set up is the ± one count 
and the accuracy of the time base reference. 

The "Frequency Heterodyne Technique", shown in figure 9, is another 
important set up used when the oscillator to be studied is more stable 
than the frequency counter time base. One must use another oscillator as 
the reference. First, the oscillator frequency is translated down to a 
value which can be conveniently measured by the counter The translation 
is realized by mixing the oscillator signal with the reference signal set 
at a convenient different frequency and by detecting only the difference 
frequency (beat frequency) . A synthesizer can be very helpful in many 
cases. The optional link between the reference and the counter allows an 
increase in the counter performance and a precise measurement of the ab- 
solute frequency. Statistical analysis is done on the beat frequency as 
described for the former set up. One must refer the fluctuations to the 
nominal frequency of the oscillator. 

If the oscillator and the reference frequencies are very close, 
one may multiply each one by a different factor, creating a more sizeable 
beat frequency. Such a system is shown in figure 10. Frequency counting 
and statistical analysis are achieved as in the previous set up. Unfor- 
tunately, noise may be introduced in the multiplier chains which set li- 
mits to this technique. 

Measurements of stability in the frequency domain require the de- 
tection of phase or frequency fluctuations. While equipped to do so, 
the same set up can be used to do measurements in the time domain. These 
ex; arimental systems are called "Phase Locked Reference Systems". A , 
voltage controllable reference oscillator is phase-locked to the oscilla- 
tor signal. When loosely locked, as indicated in figure 11, the phase 
detector delivers a signal which is proportional to the phase difference 
between the two oscillators. If the reference is considered more stable 
than the oscillator, the fluctuations of this signal is attributed to 
the random changes of the oscillator phase. When this signal is proces- 
sed through a differentiator, a low frequency signal, proportional to the 
frequency fluctuations, is obtained. It is then possible to use a volt- 
age to frequency converter driven by this signal, to generate a low fre- 
quency oscillation, fluctuating in the same manner as the frequency of 
the original oscillator. Average frequencies, over time intervals t , 
are determined by a counter and statistical analysis is performed as in 
the previous systems. , 

If the reference oscillator is tightly locked in phase to the 
oscillator studied, as shown in figure 12, the command signal, applied 
to the reference is proportional to the frequency changes between the 
two oscillators. Again if the reference oscillator is much more stable 
than the other oscillator, this command signal f 1 actuates in the same 
way as does the frequency of the oscillator. Statistical analysis of 
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this signal Is performed with the set up described in the loose phase 
locked reference case. 

The two systems act simply as frequency translators with non unity 
conversion factor. The sensitivity of these systems is enhanced by a 
factor proportional to the ratio of the nominal oscillator frequency to 
the nominal frequency of the converter. Time domain regions, where the 
results are significant, are determined by the servo loop characteristics. 
Usually a second order loop is used in order to achieve optimum perfor- 
mance; the parameters to adjust are the natural frequency and the damping 
factor [27] , [28] , [29] . The loose Phase Locked Reference System is nor- 
mally used when frequency stability for averaging times below 1 sec is 
needed, while the tight Phase Locked Reference System is preferred for 
longer averaging times. 

Obviously these systems are well suited for the measurement of phase 
or frequency fluctuations in the frequency domain. For this type of mea- 
surement, it is necessary to measure the spectral density at the phase 
detector in the first system or the spectral density of the command signal 
in the second system. Since for highly stable oscillators the information 
lies in a spectrum containing frequer es much lower than one hertz, a 
very low frequency spectrum analyser ^.s required. Digital real time spec- 
trum analysers are indicated for such measurement but they are expensive. 
The method of bandpass filtering proposed by Rutman can be applied to give 
measurements either in the time domain or in the frequency domain with 
these systems [ 5 ] . 

All the experimental systems described previously deal with the 
measurement of frequencies or periods averaged over a finite observation 
time, t . Characterization of frequency stability is done through statis- 
tical analysis on an ensemble of these averaged frequencies. I*- *as 
pointed out in the theoretical section that the frequency stab_x^ty of 
an oscillator can also be characterized by a measurement of phase diffe- 
rences. To do so. one has to measure the phase of an oscillator and to 
calculate phase differences, spaced in time by an interval, t . Measure- 
ment techniques of phase are well developed in the field of time scale 
implementation since a time scale can be graduated in terms of phase 
with 2 tt radians as a unit of time (one period). Thus, measurements of 
phase differences correspond to time difference measurements. 

We will describe two systems for measuring frequency stability by 
t^me differences. They are the Dual Mixer Time Difference System [30] 
and the Phase Modulated Phasemeter [3l] . Both systems give about the 
performance of ±1 picosecond time interval resolution and a limit of 
resolution of roughly 1 x10” 13 t - 1 when used to measure fractional fre- 
quency stability. Since these systems measure the phase of oscillators, 
the evaluation of the two-sample variance can be done without dead time 
and no correction factors are needed. Furthermore the systems can be used 
to compare oscillators of exactly the same nominal frequency (clock osci 1 
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lators) . 


The Dual Mixer Time Difference System is illustrated in figure 13. 
The two oscillators are of the same type and have the same nominal fre- 
quency, v D . A common oscillator translates down to each oscillator 
frequency in two identical channels. Phase comparison is done between 
the two beat signals. Zero crossing detection on each signal is achieved 
and serves to trigger pulse generators. The time interval between each 
pair of pulses is a measure of the relative phase of the two oscillators; 
the system acts as a sophisticated phase detector. Any fluctuations of 
this time interval can be considered as fluctuations of the phase of one 
or both oscillators. A time interval counter me aures these phase (time) 
differences and statistical analysis is done according to prescribed theo- 
retical calculations. 

The translation has the effect of increasing the system resolution 
by a factor which is approximately vq/v^ When used to characterize the 
frequency stability in the time domain, “the sampling time is the period 
or a multiple integer of the beat signal. Consequently a synthesizer 
used as common oscillator becomes a very convenient tool. Noise contri- 
bution from the common oscillator is greatly reduced when the nominal 
phase shift between the two boat signals is small. An adjustable phase 
shifter is then placed in series with one oscillator in order to satisfy 
this condition. (See the Appeudix for a calculation on that system.) 

The Phase Modulated Phasemeter System is illustrated in figure 14. 
Again two identical oscillators are compared in phase. Each signal fre- 
quency is multiplied, then mixed and filtered to get a beat signal. A 
non zero frequency beat is generated by modulating the phase of one oscil- 
lator with a low frequency signal. Zero crossing detection of this signal 
gives pulses which can be located in time when compared to a reference 
signal triggered by the modulating signal. This time interval is a mea- 
surement of the phase of one oscillator compared to the other. If one 
oscillator is considered as a reference, the time interval will be a raea- 
suxement of the phase of the other oscillator. As in the case of the 
Dual Mixer Time Difference System, the frequency stability is calculated 
from this ♦'ime interval recording for different observation times by 
simple phase differences and statistical weighing. With this system, the 
measurement resolution is increased by the up conversion factor. 


The various measurement systems described show a certain hierarchy 
in the parameters evaluated. When one has access to a signal proportional 
to the phase, he is allowed to calculate any combination of phase diffe- 
rences, then, any variances. In fact, it is possible to calculate all the 
structure functions of the phase fluctuations which gives a complete mea- 
surement of the frequency stability In the time domain. This is possible 
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with the two last systems and the loose Phase Locked Reference System 
when the phase is recorded instead of being differentiated. All the other 
systems described give access only to frequency and its fluctuations, and 
this limits the amount of statistical information which can be obtained. 
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APPENDIX: EXPRESSION OF THE FREQUENCY STABILITY IN THE TIME 

DOMAIN FOR THE DUAL MIXER TIME DIFFERENCE SYSTEM 



¥ 


* 


The phases at the output of oscillator 1 and oscillator 2 are : 
$ 1 (t) - u^t + , 

and ^(t) “ ai 2 t + 0 2 + ^2^ * 

where and are the angular frequencies, 9^ and are-, the 

initial phase offsets and <4^(t) and (^ ( t) are the phase fii - u 1 ions. 
The phase at the output or tne common oscillator is : 


$ R (t) 11 'J R t + e R + <t? R ;t) 

with the corresponding meaning for each parameter or variable. 

At the input of each mixer, the phase of the oscillator signals 
can be represented by che eq rations already given where 9 includes the 
phase shift added by the variable shifter and ^.(t) and ^(t) include 

1 i. 

any phase fluctuations added b' the transmission links. The phase of each 
reference signal can be wrxtten as : 


\,1 M • V + 9 R,1 +> R.l (t) 


V2 (t> “ “r 1 + 9 R,2 + 'R,2 U> 


in order to recount for any phase shifts and phase fluctuations introduced 
by the isolation amplifiers and the transmission links. We T '’i 1 1 ac e below 
how these phase perturbations can be made negligeable in the measurement 
system. 
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The output of each mixer deliver'' low frequency beats whose phase 


% 


< 

HK»- 


are : 



(e R,l-V 


\ 

U 

B,2^ “ 

< l ■ ,2 - ®2> 

+ V R,2^ ' 

1 


These two signals have s ime evolution represented as follow : 



of the relative phase between the two oscillators. 

At a certain time, ^ , the phase of the first beat signal is such 
that the signal is zero; then 


(“R - ^l) t l,i + 0 R,1 " 9 1 + ^R,l (t l,i ) ~ ' t ’l (t l,i ) “ m2ff 


At a certain time t ^ ^ , the phase of the second beat signal corresponds 
to the same criterion : 


c-9 
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-rtVo ' ■-« 


(w R~ W 2 )t: 2,i + 6 R,2“ e 2 + *R,2 (t 2,i* " *2 (t 2,i* “ (m + n)2lT 
where m and n are integers. 

Now if we impose that: = u>^ and define * w R -w^ , we 

obtain the identities 


m 

■-^rs 


t . 


a, B t 2,i + 0 R,2“ 6 2 + ^R f 2 (t 2,i ) '^ 2 ^ = (m + n) 2ir 


“8*1,1 + 0 R,1~ 0 1 + ^R,l (t l,i ) " ^ (t l t i* = m 2lT 


and by subtraction : 


u) B (t 2,i - t l,i ) + 0 R,2 -0 R,l~ 9 2 + e i + '^R,2 (t 2,i ) ” ^R.l^l.i* 


-^? 2 (t 2 > i) ~ = n 


2ir 


In this expression, 0 D 0 and 0 D .. are constants different by a value 
introduced by the isolation amplifiers and cable lengths. It is possible 
to define a constant phase offset. 


A0 = 0 R,2 - 0 R,1 - 0 2 + 6 1 


which can be adjusted by the phase shifter. By doing so, we also set the 
nominal value of t_ - t 1 , ; for small time offsets, the two phase flue- 
tuation terms coming from the reference oscillator are correlated and their 
difference is negligeable when compared to the phase fluctuation difference 
of the two oscillators. Within this approximation, the time difference 
becomes : 


At. = t. - t 


2,i l,i 


^2 (t 2,i ) ~ ^l (t l ,i^ _ A0 + n 2 tt 


m U) 

B B 
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If we look at the next pair of zero crossings we obtain, in a 
similar way, the time difference : 


At i+1 ~ *2,1+1 ~ t l,i+l 


= ^2 (t 2,i+l* ~ V 'i (t l,i ) _ AQ + n2ir 


CO Ul 

B B 


and for the following pair 


At i+2 *2,i+2 ~ t l,i+2 


= <t '2 (t 2,i+2 ) ^l (t l,i+2 ) _ A0 nZn 


“B 


W B W B 


The definition of the two sample variance in term of the second 
difference of phase is : 


a z (x) = -=■ 




2ttvot 


where the averaging time, x , is the time interval between and t^ + ^ 

or the time interval between two successive phase measurements of each 
oscillator. Such a linear combination of phase fluctuations can be 
obtained by grouping the time differences just derived. Then 

^2^2 i+2^ -2 ^2^ t 2 i+1^ + ^2^2 i^ 

At. 10 - 2At. ,, + At. = — -* 1 - 


'i+2 


i+1 


B 


_ ^ > l (t l,i+2 ) ~ 2 ^l (t l.i-H ) + ^. (t l.i ) 




If the two oscillators are statistically independent, we can write : 


(At i+2" 2At i+l 


+ At )" = [ < ^ (t 2,i+2 ) 2 ^2 ( t 2,i-tl ) + ^2 (t 2,i ) | 2 

+ ~ 2 ^l (t l.i+l ) + ^l (t l.i ) j' 


and in terms of the two sample variance of each oscillator : 





(4t 1+2 -24t 1+ i + 4t 1 ) -2^) T2[o|(T)+oJ(t)] 


This relation is the corner stone of the Dual Mixer Time Difference System. 
It shows that a linear combination of the time interval measurements give 
a value proportional to the two sample variance of each oscillator. When 
the two oscillators are identical, o^it) * o^(x) , and the calculated value 
is twice the value for each oscillator. If one oscillator is much more 
stable than the other oscillator, a^(x) « o|(x) , the calculated value is 
directly the value for the instable oscillator. In these calculations, we 
consider three successive pairs of time event; they are then spaced in time 
by an interval: x = (v ^) -1 * ^ eat period. It is also possible to skip 

a determined number of zero crossings between each pair of time events. In 
this case, the observation time is a multiple integer of (v „) _1 . 

O 

A similar type of calculations applies to the phase modulated phase- 
meter. 
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OSCILLATOR OSCILLATOR 



Figure 1. Block diagram of the ideal experimental set up used to measure the frequency 
of oscillator (1) over an averaging time r. The frequency of oscillator (2) is 
assumed to be free of fluctuations. 



' f 




A 

I'V 'U 

4 









<u 





• 

V 


3 

ki 

W) 

s 



<U 

3 


CN 


'o 

<u 


"kj 




o 


Vn 


<; 

V 

Nj 






K 




a 


fj 






"Is* 


o 




I 






274 


Figure 2. (a) Time sequence used in the measurement of the frequency of an oscillator. 

(b) Hypothetical result of the measurements of the frequency of oscillator (1) in the set up of figure (1). 
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Figure 3. Histogram of the hypothetical results of the 
frequency fluctuations diown in figure 2. 
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Figure 4. Hypothetical histograms obtained by repeating the measurements of frequency on N samples for various 
averaging time r. 
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(g) and (h): Modified sample variance. 
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Hadamard Variance . . 

Figure 6 (Cont.) 
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HETERODYNE TECHNIQUE 
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COMPARATOR 
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Figure 10. Block diagram illustrating the Offset Frequency Comparator method. 



PHASE LOCKED REFERENCE ( L oose ) 
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Figure 11. Block diagram illustrating the Phase Locked Reference System. In this system the reference oscillator 
is loosely locked. 
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Figure 12. Block diagram illustrating the Phase Locked Reference oscillator in the 
tightly locked mode. 
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DUAL MIXER TIME DIFFERENCE 



Figure 13. Block diagram illustrating the Dual Mixer Time Difference System. When the two 
oscillators have tendency to move off frequency they must be phase locked 
loosely together. This loop is not shown on this figure. 
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PHASE MODULATED PHASEMETER 



Figure 14. Block diagram illustrating the Phase Modulated Phasemeter System. 
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Table 1 

Asymptotic Behaviour of Various Variances for Power Law 
Spectral Models and Polynomial Drifts 


ASSYMPTOT/C BEHAVIOUR of VARIANCES 
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QUESTIONS AND ANSWERS 

DR. GERNOTM. R. WINKLER, U. S. Naval Observatory: 

I have nothing to criticize. I think it was a very clear explana- 
tion. However, I think there is still a problem if you want to in- 
troduce a generally educated person to the subject because of the 
terminology, which was, of course, adopted a long time ago; it is 
misleading. And I would like to suggest an additional approach to 
such a generally educated person. 

There are two ways, two main distinctions, in which we can 
measure and/or characterize frequency instability. The first one, 
called time domain measurement, essentially measures and/or charac- 
terizes the instabilities by measuring phase differences. And we 
obtain a statistical distribution of the carrier frequency, and we 
characterize that statistical distribution of the carrier frequency. 

In the frequency domain, we interpret the variations and meas- 
ure them as variations of the modulation frequency or of a modula- 
tion frequency around a fixed carrier. 

Now, I think that this is the first thing. Now in the time do- 
main, there are again two essentially different methods to do it. 
One is to obtain samples of the time readings which are equally 
spaced, and then you look at the statistics and have various ways to 
characterize it. 

The other way essentially is counting phase differences between 
zero crossings, and you obtain your desired statistics this way. 
Now by doing that you inevitably have the question of whether you 
have dead time or not dead time. And you have the various variations 
of your methods. 

But I believe the essential point which we tried to get across 
is that the distinction of time domain or frequency domain is not in 
frequency, but it is the distinction of frequency measurements of a 
carrier or frequency measurements of a modulation frequency around a 
fixed carrier. 

In both cases, we really measure frequency. But they mean two 
different things. Thank you. 

DR. HARRY PETERS, Sigma Tau: 

I also thought that was one of the best summaries that I have ever 
heard on the subject. If you will allow me, I would like to make 
one comment regarding the use of such systems. 

In the use of any of these systems, you need a frequency ref- 
erance. And preferably the reference should be much better than the 
things you wish to measure. If you wish to measure a crystal in a 
rubidium in 1 to 100 seconds averaging time, of course you have no 
standard which is significantly better. And for such systems, the 



use of a cesium as a reference in all of these systems is particu- 
larly inappropriate because you need a subsidiary standard as sort 
of a flywheel for all measuring times out to on the order of a thou- 
sand seconds. And this is why many people want hydrogen masers in 
their systems. For any of these systems, you would like a standard 
which is superior to all the other standards for all the measuring 
times in which you are interested. I had many other extensions of 
these comments, but I think I will stop here. Thank you. 

MR. DAVE ALLAN, National Bureau of Standards: 

Let me clarify one thing the novice to the community. I think they 
have done an outstanding job in showing how you can characterize an 
oscillator as to the random fluctuations basically. One must be 
very careful; to really characterize an oscillator, there are all 
kinds of systematic effects that must be included as far as the 
manufacturer and the user are concerned. The dependence upon 
temperature, pressure, humidity, whatever you have, that the oscil- 
lator may depend upon, is a whole set that must be included in a 
proper characterization of an oscillator. 



